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ABSTRACT 



Context. Shell-type supernova remnants (SNRs) exhibit correlations between radio surface brightness, SNR diameter, and ambient 
medium density, that between the first two quantities being the well known Y.-D relation. 
Aims. We investigate these correlations, to extract useful information about the typical evolutionary stage of radio SNRs, as well as 
to obtain insight into the origin of the relativistic electrons and magnetic fields responsible for the synchrotron emission observed in 
' radio. 

Methods. We propose a scenario, according to which the observed correlations are the combined effect of SNRs evolving in a wide 
O !■ range of ambient conditions, rather than the evolutionary track of a "typical" SNR. We then develop a parametric approach to interpret 

the statistical data, and apply it to the data sample previously published by Berkhuijsen, as well as to a sample of SNRs in the galaxy 
M33. 

Results. We find that SNRs cease to emit effectively in radio at a stage near the end of their Sedov evolution, and that models of 
synchrotron emission with constant efficiencies in particle acceleration and magnetic field amplification do not provide a close match 
to the data. We discuss the problem of the cumulative distribution in size, showing that the slope of this distribution does not relate to 
the expansion law of SNRs, as usually assumed, but only to the ambient density distribution. This solves a long-standing paradox: the 
almost linear cumulative distribution of SNRs led several authors to conclude that these SNRs are still in free expansion, which also 
implies very low ambient densities. Within this framework, we discuss the case of the starburst galaxy M82. 
' Conclusions. Statistical properties of SNR samples may be used to shed light on both the physics of electron acceleration and the 

' evolution of SNRs. More precise results could be obtained by combining data of several surveys of SNRs in nearby galaxies. 

Key words. ISM: supernova remnants - Methods: statistical - Radiation mechanisms: non-thermal - Acceleration of particles - 
Galaxies: individual: M 33, M 82 
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. 1 . Introduction 12005b have also argued that this correlation may be affected by 

^ selection effects. For Galactic SNRs, Green d2005t highlighted 

Radio emission is a quite common property of shell-type super- the selection effects; (1) in the surface brightness, with a com- 

. _h . nova remnants (SNRs). The intensity of the (synchrotron) radio pleteness limit of about 10 -20 w m -2 Hz -i sr -i > while SNRs be _ 

^ ■ emission is related to the magnetic field strength and the amount low this Hmit are predominantly in regions where the Galactic 

H i of accelerated electrons. However, the mechanisms leading to background is low; (2) in the angular size , so that young but dis . 

J*? , both the magnetic field amplification and the electron injection t£mt SNRs may be missed For SNRs in Qther galaxieSi Urogevi6 

at the SNR shock and their respective efficiencies remain poorly et al {2Q05) a i so showed that there may be a selection effect 

constrained. To investigate these processes observationally, there in the integrated flux (valid for unreS olved or mildly resolved 

have been detailed studies, mostly in X-rays, of some selected SNRs) wkh the effect of leadi tQ an observed slope shallower 

SNRs (see e.g., Cassam-ChenaT et al. 120071 Cassam-Chenai' et t ba jj t b g intrinsic one 
al. 120081 and references therein). Relevant, complementary in- 
formation should also be extracted from a statistical analysis of 

SNR data samples. Even though selection effects may affect this correlation to 
On the observational side, a well-known (even though not some level, we believe that the relation itself has a physical ori- 
widely accepted) statistical relation is the so-called "£-D re- gin, and that one can therefore extract from it information about 
lation", namely the empirical correlation discovered between the processes involved in the SNR radio emission: therefore, un- 
file SNR size (D) and its radio surfac e brigh tness (£). Various derstanding its origin could eventually contribute to constrain- 
authors (e.g., Clark & Caswell 119761 M ilne 119791 Cas well & ing the efficiency of these processes. There have been several 
Lerche |199l1 Case & Bhattacharya |1998l Urosevic et al. 12005b attempt s (e.g Shklovs ky [T9601 Van der Laan [T962l Poveda & 
have investigated this correlation. Originally applied as a tool Woltj er [T968l Kesteven[l968 or more recently Duric & Seaquist 
for estimating SNR distances, it has been found to be unreliable 1986 Berezhko & Volk 2004) to explain the £-D relation as the 
for this purpose. Some authors (e.g., Green 2005 Urosevic et al. average evolutionary track of a "typical" SNR. In all of these 
cases, the slope of the correlation is assumed to correspond to 
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Fig. 1. Distribution of SNRs from Berkhuijsen (1986) sample in 
the lg D-lg Z parameter plane (D is measured in pc, while Z is in 
Wm^Hz^'sr 1 ). The dot sizes are proportional to lgn values, 
the legend showing in order sizes corresponding to lg n Q from 
-1.5 to 1.0, in steps of 0.25 (n being measured in cm -3 ). It is 
apparent from this figure that smaller (and brighter) SNRs are 
typically located in a denser medium. 



However, there is clear evidence that this assumption is in- 
correct. Berkhuijsen ( 1986) has found tight correlations of both 
Z and D with the ambient density (n ), in the sense that smaller 
(and brighter) SNRs are typically located in a denser medium 
(indeed, the correlation between n Q and D had already been 
known for several years; see e.g., Fig. 4 in McKee & Ostriker 
119771 1. The best-fit results given by Berkhuijsen are 



D 



15n -0.39 ± 0.04 pc 



Z =* 6x 10 



-20, 1.37±0.21 



Wm^Hz-'sr 



(1) 
(2) 



Berkhuijsen concluded that the Z-D relation (Z oc D^) is just a 
secondary effect, while the two primary relations are those of D 
and Z with « D . 

Although this conclusion may appear rather extreme, it is 
quite obvious that the Z-D relation contains SNRs that evolve 
in very different ambient conditions. We then assume that the 
correlations between Z, D, and n do not directly reflect the evo- 
lution of a "typical" object, but are rather the combined effect 
of the evolution of different SNRs expanding in different ambi- 
ent conditions. That the correlations with n are statistically so 
tight also suggests that D and S are more sensitive to the ambi- 
ent conditions than to the SNR evolution. The combined effect 
of these correlations is that a clear trend of n across the Z-D 
relation is observed (see Fig. [TJ, namely that smaller SNRs are 
preferentially located in higher-density environments. 

In principle, correlations of Z and D with n , accounting for 
at least qualitatively the trend of n a across Z-D relation, could 
also appear in the case of the evolution of an individual object, 
provided that it expands in a medium with a highly structured, 
fractal, density distribution. This would reflect that, during its 
life, a SNR always preferentially expands towards the direction 
in which the ambient density is lower. Therefore, at any time, 
the "effective" ambient density would be close to the "lowest" 
ambient density in the volume occupied by the SNR. In this way, 
one could explain why large SNRs typically seem to expand in 
a low-density medium. 

However, it is hard to justify, in this scenario, the absence 
of small SNRs in low-density media (which would be the case, 
when the supernova itself is located in a low density region), as 
well as that low ambient density values are measured for all ex- 



tended SNRs (which requires that such low density regions are 
ubiquitous in the Galaxy, on scales of tens of parsecs or even 
less). In addition, if the effects of the fractal interstellar medium 
were dominant, virtually all SNRs should have a much brighter 
limb on one side, which is not observed. Therefore, for all these 
reasons, we conclude that the "fractal ambient density" hypoth- 
esis is implausible, and we do not consider it any longer in this 
paper. 

One of our goals is to show that the best-fit line usually re- 
ferred to as the "Z-D relation" provides only a minor part of 
the information present in the data, while additional informa- 
tion could be extracted by analyzing in detail the distribution of 
points in the Z-D-n parameter space. For this reason we pro- 
pose a rather general (parametric) scenario, with the aim of con- 
straining the physics of the electron injection and magnetic field 
behavior in SNR shocks and/or the SNR evolutionary phase in 
which they are most likely to be observed in radio. 

The plan of the paper is the following. In Sect. 2, we present 
the basic ideas, assumptions and formulae that will be used in 
the rest of the paper. Section 3 is devoted to a re-analysis of the 
data of Berkhuijsen ( 1986), using the criteria introduced in Sect. 
2, and extracting constraints on the dependence of the particle 
and magnetic field efficiencies on the basic shock parameters. 
For the sake of comparison, we also apply the same technique to 
a sample of SNRs in galaxy M 33. Section 4 shows how the cu- 
mulative distribution of SNRs with diameter may be interpreted 
within our framework, and how we can explain the puzzling lin- 
ear trend of this distribution without requiring, as usually done, 
that SNRs expand linearly to large sizes. Section 5 presents our 
conclusions. 



2. Basic ideas, assumptions, and formulae 

Our analysis is based on the fundamental criterion that the ob- 
served correlation in the Z-D parameter plane originates from 
the combined effect of evolutionary tracks in very different ambi- 
ent conditions. In this section, we implement this idea by adding 
some derived / secondary assumptions that will allow us to de- 
velop a more general scenario, on which our subsequent statisti- 
cal analyses will be based. 



2. 1 . When radio SNRs are preferentially seen 

A preliminary consideration is that the conditions in which a 
given object is most likely observed, during its evolution, are 
those in which it spends most of its time. Since SNR expan- 
sion decelerates during most of their lifetime, it is statistically 
more likely to find them when their size is close to its final 
value. We are interested in finding SNRs that are visible in ra- 
dio. Therefore, it is more important, in this case, to determine 
the evolutionary stage at which the processes responsible for en- 
hancing magnetic fields and/or for producing high energy elec- 
trons are no longer efficient. In the following, we refer to this 
phase as the "final stage" of a radio SNR, but it should be clear 
that it is not the maximum size that a SNR can reach dynami- 
cally, before merging into the ambient medium. 

We parametrize the SNR expansion by a power law (D oc 
t '°; a > 1) up to a maximum size (D2) beyond which the SNR 
is no longer detectable in radio. During its evolution of a given 
SNR as a radio source, the probability of being observed with a 
given size D is proportional to dt/dD, namely 



P(D) = aD a - l /D" where D < D 2 , 



(3) 



R. Bandiera & O. Petruk: A statistical approach to radio emission from shell-type SNRs 



3 



(the initial diameter of this evolutionary phase, D\ , not being rel- 
evant provided that (Di/Z^) 0-1 <k 1), so that the average value 
and standard deviation of the (decimal) logarithm of D are 



<lg£>> 

CT\ g D 



lgZV 
l 

fllnlO' 



1 



alnlO' 



(4) 
(5) 



For instance, during the adiabatic (Sedov) phase a = 5/2, so 
that (lgD) = lgD 2 - 0.17 and cr lgD = 0.17, while, in the later 
radiative (pressure-driven snowplow) phase, a = 7/2, so that 
(lgD) = \gD2- 0.12 and <x leD = 0.12. This means that, on 
average, SNR diameters should be rather close to D2, and that 
their dispersion should be rather small, i.e., an individual SNR, 
during its evolution, is seen to migrate only slightly in the 2- 
D parameter plane. For this reason, we propose that selection 
effects, while important to determining the overall distribution of 
points across the X-D plane, should only have a marginal effect 
on the observed probability P(£>). 



2.2. The end of the radio phase 

We consider, in particular, the end of the Sedov phase. 
According to Truelove & McKee (119991 ), it should correspond 
to a size 



D B ~28 (£ S N/10 51 erg) 2/7 «- 3/7 pc, 



(6) 



(where £sn is me energy of the supernova explosion), while the 
dynamical end of the SNR, namely where it merges with the 
ambient medium, can be placed at much larger sizes. The above 
formula has been obtained by approximating the plasma cool- 
ing function with a power law A oc r -1 ' 2 , where T is the gas 
temperature (while other papers use different power-law approx- 
imations; for instance, Blondin et al. 1998 use A oc r -1 ). For a 
generic A oc T~° relation one may find that 



n p (3+2a)/Ql+6or) -(5+2a)/(ll+6or) 



(7) 



For a changing from 1/2 to 1, the exponent of n in the above 
formula changes from 0.429 to 0.412: namely, the numerical 
value of that exponent is very weakly dependent on the power- 
law approximation used. In the following, we shall then use, 
without loss of generality, the formula (here Eq. [6} of Truelove 
& McKee ( [19991 . 

The correlation found by Berkhuijsen between D and n Q (Eq. 
Q]) is consistent with 



D-0.54(£ SN /10 51 erg)~ 2/7 D 



(8) 



namely with a D/Dg ratio that is independent of n . 

This can be also seen from Fig. [2] where the data points from 
Berkhuijsen (1986) are displayed together with Dg, as evaluated 
for £sn = 10 51 erg (solid line): we note that the line is not a fit, 
namely there are no free parameters to tune. This indicates that 
most of the known radio SNRs are observed close to the end of 
the Sedov phase, and that in general SNRs must extinguish their 
radio emission somewhere close to the end of their Sedov phase. 
It is then reasonable to use a Sedov law (a = 5/2) to approxi- 
mate the expansion law during the final phases of radio SNRs: 
therefore, in the following, whenever a numerical value for a is 
required, we shall use the Sedov value. For a Sedov expansion, 
the dashed line in Fig. [2] applies to the value of (lg D), given by 




Fig. 2. Distribution of SNRs (from Berkhuijsen d 19861 ) sample) 
in the lg n -lg D parameter plane. For comparison, the theoret- 
ical line corresponding to the end of the Sedov phase (as from 
Truelove & McKee 1999 ) is shown (solid line), as well as that 
of (lg£>), as from Eq.|4](dashed line). 

Eq. E] The best-fit level (Eq. [8]l is only 20% lower, providing 
a good argument for SNRs being (statistically) mostly visible 
around the end of their adiabatic stage. 

It remains to be understood for which physical reason the end 
of the Sedov phase should roughly correspond to the switching 
off of the radio emission. There is also evidence that the extin- 
guishing transient evolution must be rather rapid. Otherwise, we 
should also see SNRs with X values considerably lower than that 
derived from the "L-D correlation for the same size; namely, in 
the YZ-D parameter plane, we should have points spread over the 
half-plane below the main correlation (of course, limited to the 
region of the 1Z-D plane for which one expects SNRs to be de- 
tectable). This latter piece of observational evidence is also not 
easy to explain. The underlying problem is that the physical pro- 
cesses behind the injection of electrons are poorly understood 
in the general case, and are probably even harder to model for 
conditions near to marginal efficiency. 

2.3. The "final-stage" approximation 

We introduce in a parametric form a basic set of equations to 
describe the observed correlations and distributions, and eventu- 
ally provide some constraints on future physical models of the 
injection of electrons in SNRs. That the Y.-D empirical relation 
is a power law suggests (even though it does not strictly imply) 
that all formulae of interest can be approximated by power laws, 
thus simplifying considerably the treatment. 

As a starting point, we consider the extreme approximation 
that each of them is observed very close to its final stage as a 
radio SNR, namely that each individual evolutionary track in the 
2-D plane can be assimilated to just one point, corresponding to 
its final position (D2, £2). We also assume that the dependence of 
both quantities on the ambient density (n Q ) can be approximated 
by the following power laws: 



D 2 {n ) = Kyn'S, 
2 2 (« ) = K 2 n n . 

In this limiting case, the slope of the "L-D relation 
E 2 (£>2) = K 3 Di 



(9) 
(10) 

(ID 



would simply be £ = n/m. In the above formulae, the func- 
tional dependence on other physical parameters is not given 
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explicitely; however, other quantities may be involved. For in- 
stance, if SNRs really are efficient radio emitters only until the 
end of the Sedov phase (as suggested by Fig. 2), D 2 should also 
depend on the SN energy (see Eq. [7j. 

To describe the distribution of points along the correlation, 
one must also introduce the function ft(n ), giving the probabil- 
ity of finding a SNR in a region of a given density: this probabil- 
ity combines the density distribution of the interstellar medium, 
the dependence of the SN rate on the local density, and how the 
lifetime of a radio SNR depends on the ambient conditions. For 
the sake of simplicity, and in the absence of any observational 
evidence against it, we also approximate this function by a power 
law, namely 

P(n ) = £ 4 <. (12) 
This distribution is used in Sect. 15.11 

2.4. Introducing the SNR evolution 

From this point on, we shall remove the "final-stage" approxi- 
mation, introduced in the previous section. Nonetheless the evo- 
lution of individual SNRs will still be treated in a very simpli- 
fied way, by assuming that SNR evolution in different ambient 
conditions differ only by a scaling law and, in practice, by only 
adopting power-law behaviors. 

We introduce the SNR expansion in the following parametric 
form: 

t(D,n ) = K 5 D a n b (13) 

(for instance, in the Sedov case a — 5/2, b — 1/2, and K5 oc 
E s ^ 2 ). As for the evolution in surface brightness, by assum- 
ing that a scaling law holds, the surface brightness could be ex- 
pressed in a rather general form as 

E(D, no) = /(D/D 2 (n ))S 2 (« ), (14) 

where fix) must vanish at x > 1 . In this paper, we use a power- 
law approximation 

£(D,n ) = (— ^— Y Ie(n ) = KeLTnl for D < D 2 , (15) 
\D 2 (n )l 

where = K 2 jK p { and q = n—mp. Parameters p and q in Eq.[T31 
can be derived independently. If the data sample is not heavily 
affected by selection effects, these parameters can be evaluated 
by simply applying a bilinear regression. 

Determining the functional dependence of E(Z), n„) would 
also allow one to derive the trajectory of individual SNRs in the 
parameter plane. They are simply given by the function £(£>) 
for a constant value of n , which in principle differs from the 
2(D) relation, as traditionally obtained, because the latter rela- 
tion is obtained by combining cases with different n values. In 
the power-law case, the slope of the evolutionary track of an in- 
dividual SNR is then given by the exponent p. In the following, 
we shall provide some evidence that the value of p is different 
from that of f: this means that evolutionary tracks in the ~Z-D 
parameter plane have a different slope from that of the overall 
Y.-D relation. 

Also m can be derived, by fitting Eq.|9]to the D-n Q data (pro- 
vided that (D) is a constant fraction of D 2 , as from Eq.|4j. On the 
other hand, there is no way of deriving a and b (defined in Eq.fTSl 
directly from the correlations between 2, D, and n . Coefficient 
a could be inferred, in principle, only by studying the distribu- 
tion of points about the main correlation (a point that will be 



treated in a forthcoming paper), while there is no way of esti- 
matig the exponent b, because it would only affect the distribu- 
tion of points with n along the correlations, a piece of informa- 
tion that is already included in the definition of the distribution 
P (Eq.CTJ). 

In principle, a distribution of SN energies may also con- 
tribute; but, for the present analysis, we assume that SN energies 
are not correlated with any other quantity and therefore that a 
distribution of energies would just produce an additional disper- 
sion across the correlations, without affecting any of the above 
slopes. Some correlation could be possible, in principle, if dif- 
ferent population stars have different distributions of their SN 
energies. However, to our knowledge no evidence in favour of 
this has been presented so far. 

3. Statistical analysis of the data 

3.1. Data sample and best-fit parameters 

We now apply the analysis outlined above to the data published 
by Berkhuijsen (I1986l l. That paper presents a fundamental work 
on the subject and, although since then a great number of sur- 
veys of higher accuracy have been performed, it still contains 
the most extended data sample of SNRs in which, in addition to 
the SNR radio surface brightness and size, quantities derived by 
other spectral bands are also tabulated. We extracted from this 
data sample all SNRs with available data on E, D (radio), and 
n as well. We excluded SN 1006, because it is now known that 
the bulk of its X-ray emission is non-thermal. The total num- 
ber of selected objects is 34: the original data sample is given 
in Table 1 (for two objects, Cas A and Tycho, in which two dif- 
ferent values of n Q are given, we took their geometrical mean). 
Since distance estimates for Galactic SNRs have changed with 
time, in the last column of Table 1 we list the SNR sizes obtained 
from the most recent version of the Galactic SNR catalog, by 
Green ( 120091 ): they differ substantially from Berkhuijsen's val- 
ues only for Kepler and Vela. The estimated distances of the 
Large Magellanic Cloud (LMC) and the Small Magellanic Cloud 
(SMC) have also slightly changed, from 55 kpc and 63 kpc (as 
in Berkhuijsen 1986), to 48 kpc and 61 kpc (Macri et al.[2006 
Hilditch et al. 2005 ), respectively. In our calculations, we used 
all of these new distances, and we revised accordingly the SNR 
linear size and density estimates (being n Q oc dr 1 ^ 1 ). The last two 
columns of Table 1 provide the D and n values that we used. 

The average ambient density is estimated by Berkhuijsen 
(1986) in a simple way, using the following relation (derived 
from LongH983]): 

n = {6ln?l 2 e- l ' 2 fl 2 lJfD-V\ (16) 

where e (taken to be 3 x 10~ 23 erg crrT 3 s _1 ) is the specific emis- 
sivity, / is the filling factor (taken to be close to unity), and Lx 
is the X-ray luminosity. While the exact values of e and / are 
not important to our statistical analysis, it is crucial that these 
quantities remain constant, or at least independent of other pa- 
rameters (such as size and surface brightness). In spite of its 
simplicity, this formula provides reasonably good results. Upper 
limits to the uncertainty in this n Q evaluation can be derived from 
the dispersion about the lgD-lgn regression. Based on the as- 
sumption that the measured dispersion depends only on the un- 
certainties in n Q , one obtains <x(lgn ) = 0.42; while also tak- 
ing into account the dispersion in D, as modeled by Eq. [5] for 
the case of Sedov expansion, one derives a residual dispersion 
<r(lg« ) = 0.18, namely a typical uncertainty in the density de- 
rived by Berkhuijsen of only about 50%. 
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Table 1. The data sample (extracted from Berkhuijsen 1986) 





Object 


-lgS 


D 




D 










(as from 


(after distance 






(Wrrr 2 


Berkhuijsen) 


revision) 






Hz-'sr 1 


) (pc) 


(cm" 3 ) 


(pc) 


(cm" 3 ) 


GAL 


W44 


19.409 


26.0 


0.22 


25.0 


0.22 


GAL 


Cas A 


16.721 


4.1 


4.94 


4.9 


4.52 


GAL 


Tycho 


18.921 


5.4 


3.64 


5.6 


3.57 


GAL 


RCW103 


19.377 


9.6 


5.04 


9.0 


5.21 


GAL 


Kepler 


18.481 


3.8 


7.70 


2.5 


9.49 


GAL 


W49B 


18.432 


11.0 


0.70 


10.1 


0.73 


GAL 


VelaXYZ 


20.367 


36.0 


0.02 


18.5 


0.03 


GAL 


RCW86 


20.276 


35.0 


0.08 


28.1 


0.09 


LMC 


0453-685 


19.730 


20.0 


0.44 


17.5 


0.47 


LMC 


0454-665 


19.629 


16.0 


0.44 


14.0 


0.47 


LMC 


0455-687 


20.299 


51.0 


0.06 


44.5 


0.06 


LMC 


0500-702 


20.590 


30.0 


0.12 


26.2 


0.13 


LMC 


0505-679 


20.374 


9.0 


1.34 


7.9 


1.43 


LMC 


0506-680 


18.976 


12.0 


2.14 


10.5 


2.29 


LMC 


0509-675 


19.231 


7.0 


4.64 


6.1 


4.97 


LMC 


0519-697 


19.524 


24.0 


0.10 


20.9 


0.11 


LMC 


0519-690 


18.984 


7.0 


6.76 


6.1 


7.24 


LMC 


0520-694 


20.435 


37.0 


0.20 


32.3 


0.21 


LMC 


0525-660 


19.480 


30.0 


0.60 


26.2 


0.64 


LMC 


0525-696 


18.293 


22.0 


4.12 


19.2 


4.41 


LMC 


0525-661 


18.544 


16.0 


2.68 


14.0 


2.87 


LMC 


0527-658 


20.807 


61.0 


0.12 


53.2 


0.13 


LMC 


0528-692 


20.168 


30.0 


0.16 


26.2 


0.17 


LMC 


0532-710 


19.947 


43.0 


0.12 


37.5 


0.13 


LMC 


0534-699 


20.356 


33.0 


0.44 


28.8 


0.47 


LMC 


0534-705 


20.481 


37.0 


0.08 


32.3 


0.09 


LMC 


0535-660 


18.287 


12.0 


4.30 


10.5 


4.60 


LMC 


0536-706 


20.481 


30.0 


0.10 


26.2 


0.11 


LMC 


0543-689 


20.758 


45.0 


0.10 


39.3 


0.11 


LMC 


0547-697 


19.736 


50.0 


0.10 


43.6 


0.11 


LMC 


0548-704 


20.151 


22.0 


0.34 


19.2 


0.36 


SMC 


0045-734 


19.802 


27.0 


0.10 


26.1 


0.10 


SMC 


0102-722 


18.570 


8.0 


8.28 


7.7 


8.41 


SMC 


0103-726 


20.807 


54.0 


0.26 


52.3 


0.26 


By performing 


linear rej 


ressions 


between the lo| 


;arithmic 



quantities (which is equivalent to assuming constant relative er- 
rors in the measurements), we obtain: 



m = -0.37 + 0.04, 
f = -2.06 + 0.34, 
p = -0.89 + 0.57, 
q = +0.62 + 0.25, 



(17) 
(18) 
(19) 
(20) 



(where l-cr uncertainties are indicated). It is apparent that, while 
m and £ are found to have reasonably small uncertainties, the 
uncertainties of p and q are larger. The reason is that, in the data, 
there is a near-degeneracy between p and q, as is well shown by 
a plot of the confidence levels (Fig. [3). In this sense, a combined 
quantity that can be far more reliably determined is: 



4 -0.39/;- = 0.97 + 0.14. 



(21) 



A potential problem of this sample, and of SNR samples in 
general, is the presence of selection effects. In the introduction, 
we mentioned the analyses of Green (2005) and Urosevic et al. 
(2005) on this subject. The points raised by Green (2005) are 
more appropriate to our sample, which consists only of SNRs 
located in our Galaxy and the Magellanic Clouds. 

Even though a detailed treatment of the selection effects is 
beyond the scope of this paper (but will be approached in a 
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Fig. 3. Plot of the confidence levels in the p-q parameter plane. 
The levels plotted correspond to 1,2, and 3-cr confidence lev- 
els, while the line indicates the maximum spread direction (see 
Eq.[2TTi. For comparison, two theoretical predictions are plotted: 
"Model A" (-17/4, 0) indicates the case with constant efficiency 
in both particle acceleration and magnetic field compression plus 
amplification (Berezhko & Volk l20D4b : while "Model B" (-2, 0) 
refers to the case in which particles are accelerated with constant 
efficiency but the magnetic field is constant (see text). 
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Fig. 4. Plot of the l-cr and 3-cr confidence levels in the p-q 
parameter plane, for the whole sample (solid lines), as well as 
for selections of the 25 and 30 SNRs with the highest surface 
brightness (short-dashed and long-dashed lines, respectively). 



forthcoming paper), for the sake of illustration we repeated the 
computations that led to Fig. [3] but on subsamples containing, 
respectively, the 25 and 30 SNRs with the highest radio surface 
brightness with corresponding thresholds in lg E of -20.356 and 
-20.481 respectively (see Table 1). Figure|4]shows a comparison 
of the confidence levels (only l-cr and 3-cr, for figure clarity) for 
the two subsamples, superimposed on those for the entire sam- 
ple (as in Fig. [3]). It is apparent that any additional selection in 
surface brightness does not generate any substantial difference 
in the results, apart from a slight terms of broadening compati- 
ble with the lower sample size being analyzed by statistics and a 
slight shift toward lower p and q values (the latter point will be 
discussed in the next section). 



3.2. Testing the "constant efficiencies" model 

Among the various theoretical attempts to model the radio emis- 
sion from SNRs, one of the most recent and popular is that by 
Berezhko & Volk (2004). This paper assumes that the kinetic 
energy density entering the shock (mn V s 2 h , where m is the mean 
atomic mass) is converted with constant efficiencies into the en- 
ergy densities of magnetic field and accelerated electrons (e^ 
and £cr, respectively): this means, for instance, that the effec- 
tive magnetic field in the synchrotron emitting region decreases 
with time, simply because the SNR shock is slowing down. For 
synchrotron emission with a power index -0.5 (namely the av- 
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Fig. 5. Same as Fig. [3] in the g-h parameter plane. 



erage index for radio SNRs), the surface brightness should scale 
as 



E cc KB 3/2 D oc (e CR 4 l4 )D <* (« D V s 2 h ) 7/4 D 



(22) 



A further assumption of this model is Sedov expansion, which 
implies that mn V 2 , ~ Es^/D 3 , so that one finally obtains 



ZccE 7 ^D- 17 >\ 



(23) 



i.e., with p = -4.25, and q — (labeled as "Model A" in Fig. [3]). 
Namely, according to this model, individual SNR tracks in the 
E-D plane must be rather steep and "independent of the ambient 
density". Since Berezhko & Volk (2004) state that the slope of 
the X-D relation should represent the slope of individual evolu- 
tionary tracks, they predict that £ = -4.25 should be the slope 
of the E-D relation. However, neither their predicted value for 
£ matches the data, nor does (and at an even higher significance 
level) their predicted (p, q) pair (see Fig. [3). It is unlikely that 
this mismatch is a mere consequence of a sample incomplete- 
ness in surface brightness. Figure [4] shows that for subsamples 
in which a further selection in surface brightness has been ap- 
plied the barycenter of the confidence levels moves only mildly. 

A more appropriate model (in the sense that it is "only" about 
2-cr away from the best-fit model values) would be one in which 
electrons are accelerated with constant efficiency (ecu oc n V 2 h ) 
but the magnetic field is taken to be constant, not only during the 
evolution of an individual SNR but also among different SNRs. 
This happens, for instance, if the post-shock field has been com- 
pressed only by the shock, i.e., is proportional to the ambient 
field, which in turn is roughly constant (see e.g., Crutcher et al. 
2003), in near equipartition with the interstellar thermal pressure 
of the diffuse interstellar medium. This case does not exclude the 
presence of an extra field amplification, provided that it yields a 
constant factor. As in the previous case, but assuming B to be 
constant, one can now write 

£ oc KD oc <e cr D oc (n V 2 h ) D oc E SN D~ 2 , (24) 
namely p — -2 and q = (which is labeled "Model B" in Fig. [3]). 

3.3. The results with a more "physical" flavour 

We can approach the problem from the opposite direction, by 
trying to translate the information derived in terms of (p, q) into 
constraints on the physics that controls the magnetic and cosmic- 
ray efficiencies. We assume that 



where g and h are free parameters. This is not the most general 
case, but it simply relates the efficiencies to primary local quanti- 
ties encountered by the shock. Based on this assumption, and an 
expansion law t oc D" passing through the endpoint of the Sedov 
phase (Db, fg), the equation for the surface brightness becomes 

Z oc 4 4 N fl - 3)/14 D 1 -( fl - I > A «r (3fl - 4)/,/7 . (26) 
In the case of Sedov expansion, this equation simplifies into 



ZocE^D 



l/2 n l-3A/2 M g-A/2 



(27) 



Figure [5] shows the confidence levels in this new pair of parame- 
ters. An advantage is that the direction of maximum dispersion is 
almost parallel to the h axis, which means that at least the best-fit 
value for g is well determined. We have 



g = 1.25+0.14, 
h = 1.26 + 0.38. 



(28) 
(29) 



(25) 



For comparison, the "constant efficiencies" model prescribes 
that g — 1.75 and h = 3.5, as can easily be derived from Eqs.l23l 
and|27] What we have found here is not different from the previ- 
ous section, but is simply displayed in a more physical way. 

To summarize, in most SNRs the constant efficiency assump- 
tion (namely for both field amplification and particle accelera- 
tion) does not hold. This may not be surprising, in the view that 
the statistically most common cases are those of SNRs close to 
their radio endpoint, namely when particle acceleration is close 
to being halted. Using observations to test these critical cases 
may indeed be important to obtaining a clearer understanding 
of the physical processes responsible for magnetic amplification 
and particle acceleration. 

4. Results from an independent sample: M 33 

As already mentioned in the introduction, SNR samples studied 
in other galaxies provide promising input to this kind of anal- 
ysis. On the other hand, the presence of selection effects could 
affect the results (Urosevic et al. 2005), and a statistical anal- 
ysis of these samples should include a careful study of these. 
Extragalactic SNR samples will be considered in more detail in 
a forthcoming paper. For mere comparison with what has already 
been obtained using the sample of Berkhuijsen, we present here 
the results for a M 33 SNR data sample, which is a sample com- 
pletely independent from that used so far. 

The sample was obtained by selecting SNRs for which ra- 
dio fluxes are given by Gordon et al. (1999), X-ray fluxes by 
Pietsch et al. (2004), and (optical) diameters by Gordon et al. 
(fT998l . The X-ray survey of Plucinsky et al. ( [20081 , based on 
Chandra data, was also used to solve some cases of uncertain 
identification. In this way, we selected 22 SNRs: the data sample 
is shown in Table [2] The first 3 columns report the SNR iden- 
tification numbers in the various catalogs, respectively, Gordon 
et al. JT9981 labeled by "opt" in the T able) Gordon et al. ( [19991 
labeled by "rad"), and Pietsch et al. ( 120041 labeled by "xray"), 
while the next 3 columns show, respectively, the published linear 
(optical) diameters, 20 cm fluxes, and measured (i.e., absorbed) 
0.2^1.5 keV fluxes. 

To derive 1 GHz radio fluxes, we extrapolated the 20 cm 
fluxes tabulated by Gordon et al. (19993, by assuming a spec- 
tral index of -0.5. For most SNRs, Gordon et al. 0999) also 
estimate spectral indices, but the uncertainty in these estimates 
is rather large, and we therefore preferred to adopt a "standard" 
value for the spectral index. We evaluated the unabsorbed X-ray 
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Table 2. M 33 data sample 



ID 

opt 



ID 
rad 



ID 

xray 



D 

(pc) 



S(20 cm) 
(mJy) 



Flux(0.2^1.5) 

(ergcm~ 2 s _1 ) 



9 


11 


93 


18 


0.7 


4.29 x lO" 15 


11 


13 


98 


17 


0.6 


3.81 x 10~ 15 


15 


20 


106 


27 


0.6 


2.79 x 10~ 15 


20 


25 


120 


10 


0.8 


1.03 x 10~ 14 


21 


29 


121 


28 


0.9 


1.41 x 10- 13 


25 


42 


144 


27 


1.4 


3.63 x 10- 15 


27 


47 


153 


23 


1.2 


2.99 x 10 15 


28 


50 


158 


11 


0.8 


2.45 x 10~ 14 


29 


52 


161 


20 


0.5 


1.73 x 10~ 14 


31 


57 


164 


39 


1.8 


3.88 x 10~ 14 


35 


64 


179 


32 


3.5 


1.01 x 10~ 14 


42 


75 


194 


29 


0.5 


1.39 x 10~ 14 


47 


90 


207 


36 


0.2 


7.05 x 10~ 15 


53 


110 


213 


40 


0.2 


1.53 x 10~ 15 


54 


111 


214 


16 


1.3 


2.08 x 10~ 15 


55 


112 


215 


18 


4.4 


2.84 x 10~ 14 


57 


114 


220 


21 


0.4 


1.87 x 10~ 15 


59 


121 


224 


16 


0.3 


3.90 x 10- 15 


62 


125 


225 


29 


0.4 


4.87 x 10" 15 


64 


130 


230 


27 


0.5 


3.65 x 10- 15 


73 


148 


250 


17 


0.5 


1.41 x 10~ 14 


97 


181 


314 


35 


0.8 


4.84 x 10~ 15 



fluxes by taking a column density Nh = 1.0 x 10 21 cirT 2 towards 
M 33 (Plucinsky et al.[2008 ), and assuming for the average SNR 
spectrum a Raymond-Smith model with kT = 0.3 keV. Using 
WebPIMMSQ a correction factor of ~ 1 .94 is evaluated: the pre- 
cise value of this factor is not very important for our purposes, 
provided that the X-ray SNR spectra are not too different among 
themselves. Finally, a M 33 distance of 817 kpc (Freedman et al. 
1200 11 1 is assumed here; since Gordon et al. ( 1998) use a distance 
of 840 kpc, for consistency we applied a small correction to their 
published SNR sizes. 

Here are our results, to be compared with those presented 
above. The formulae equivalent to Eqs. [T7if2"0l |28l and [29] are 
respectively 



m = -0.34 + 0.07, 

£ = -2.20 + 0.46, 

p = -1.37 + 0.64, 

q = +0.52 + 0.30, 

q - 0.37/? = 1.04 + 0.20 

g = 1.31 + 0.20 

h = 1.58 + 0.42; 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 
(36) 



while Figs. 6, 7, 8, and 9 correspond, for M 33, to Figs. 1, 2, 
3, and 5, respectively. It is apparent that all results from this 
further sample show close agreement, within the quoted errors, 
with what we found above using the data of Berkhuijsen. A com- 
parison of the regression results for the two samples is given in 
Table 3. 

5. The SNR cumulative distribution with size 

5. 1 . The original paradox and how it can be solved 

The cumulative distribution of the number of SNRs with sizes 
smaller than a given diameter (N-D relation) is another statisti- 
cal distribution that has traditionally been studied. 



■H 

W 




20.5 



1.2 1.3 1.4 
lg D 

Fig. 6. Distribution of the M 33 SNRs sample in the lgD-lgl 
parameter plane (to be compared with Fig. 1). The dot sizes are 
proportional to lg n values, the legend showing in order sizes 
corresponding to lg n Q from -1.0 to 0.5, in steps of 0.25. 
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Fig. 7. Distribution of M 33 SNRs in the lgn Q -lgD parameter 
plane (to be compared with Fig. 2). 
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Fig. 8. Plot, for M 33 SNRs, of the confidence levels in the p-q 
parameter plane (to be compared with Fig. 3). 

Table 3. Synoptic table of the results of the regression analyses 



Formula 




Data Samples 






Berkhuijsen 


M 33 


D 2 (n ) oc < 


m = 


-0.37 ± 0.04 


-0.34 ± 0.07 


£2(02) oc D{ 


? = 


-2.06 ± 0.34 


-2.20 ± 0.46 


£(£>, n ) cc D"n 9 


P = 


-0.89 + 0.57 


-1.37 + 0.64 




1 = 


+0.62 ± 0.25 


+0.52 ± 0.30 


q-Xp = ii 


A = 


+0.39 


+0.37 




V = 


+0.97 + 0.14 


+ 1.04 + 0.20 


KB 312 cc n*y* 


8 = 


+ 1.25 + 0.14 


+ 1.31 + 0.20 


h = 


+ 1.26 + 0.38 


+ 1.58 + 0.42 



http : //heasarc . gsf c . nasa . gov/Tools/w3pimms . html 
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Fig. 9. Plot, for M 33 SNRs, of the confidence levels in the g-h 
parameter plane (to be compared with Fig. 5). 




1.8 



Fig. 10. Cumulative distribution, for Berkhuijsen's sample, 
where the dot sizes are proportional to the lg n value of the last 
SNR entering in the cumulative. The legenda uses for the dot 
sizes the same convention as in Fig.Q] 



For the Magellanic Clouds, Mills at al. (1979 ) derived an al- 
most linear relation (N oc D 12 ) up to sizes as large as 40 pc. 
They also argued that a similar relation should be present in 
our Galaxy. A linear cumulative distribution is usually taken as 
evidence that these SNRs are still in free expansion. However, 
when the SNR diameter is 40 pc, the swept-up mass is ~ 
(1000/n o )M o . Therefore, except for cases of exceptionally low 
ambient density, at those sizes SNRs should already be in the 
Sedov phase, and therefore strongly decelerated. A way to solve 
this paradox was suggested by Green (1984), which invokes a 
major effect of sample incompleteness. Indeed, any sample in- 
completeness is expected to substantially affect the cumulative 
distribution. Much less pronounced effects on the "L-D corre- 
lation, as well as on the correlations discussed in Sect. [3] are 
expected, since the main effect of the sample incompleteness 
should be to change the distribution of points "along" the cor- 
relations. The relevance of data sample incompleteness to the 
correlations will be discussed in a forthcoming paper. 

Within our framework, the cumulative distribution with size 
is independent of the expansion law of individual SNRs, but 
is related instead to the statistical distribution of the ambient 
medium density, as defined by Eq. Q~2] Indeed 



N(D) oc <P(n )n Q oc D' 



(l+M')/m 



(37) 



In this sense, it is not even necessary to account for the sample 
incompleteness, as done by Green ( 1984J. 

For the sake of illustration, Fig.[lO]shows the N(D) cumula- 
tive distribution of Berkhuijsen's sample. Here the close to lin- 



earity of the distribution is coincidental, since there is no reason 
to expect the sample of Berkhuijsen to be complete. Instead, the 
trend of dot sizes with SNR diameter (which is to some extent a 
different way of displaying the information contained in Fig. |2]i 
clearly shows how, for increasing size, the cumulative is more 
and more populated by SNRs located in more tenuous ambient 
media. 



5.2. The case of M 82 

Samples of SNRs in nearby galaxies have become increasingly 
available, and the close to linearity of the N—D relation is a rather 
standard property of these samples. 

A particularly interesting case is that of M 82. In this nearby 
starburst galaxy, a number of radio sources have been detected 
(Rronberg et al. 1985), which may be SNRs that are much 
brighter and much smaller in size (a few parsec at most) than 
usual. Their positions on the E— D parameter plane are in all cases 
consistent with the extrapolation to lower sizes of the best-fit E— 
D relation for other galaxies. Chevalier & Fransson (2001) pro- 
posed that these sources are SNRs expanding in a high-density 
ambient medium (with densities of order of 10 3 cm -3 ). 

Rronberg et al. (2000) also placed upper limits on the flux 
density variations in most of the radio sources: approximately 
75% of these objects are very stable, with a lower limit of ~ 
10 3 yr to their characteristic radio-emitting lifetimes. Based on 
this upper limit, Seaquist & Stankovic (2007 ) suggested that they 
may not be SNRs, but rather cluster wind-driven bubbles. Their 
main argument is that the lack of observed time variability is 
inconsistent with the estimated ages of these objects. That is, 
if they are SNRs in free expansion (with typical velocities of 
~ 10, 000 km s" 1 ), their ages should be a few hundred years at 
most, while to account for the lack of variability, velocities no 
greater than ~ 500kms _1 are required (Chevalier & Fransson 

HEED). 

A crucial point in this reasoning is the expansion regime of 
these objects, if they are indeed SNRs. A linear expansion is 
argued by both Muxlow et al. d 19941 ) and Fenech et al. (2008) 
based on the cumulative distribution with size being almost lin- 
ear. However, as we have explained above, a linear cumulative 
distribution does not imply a linear expansion, if it is caused by 
the combination of SNRs expanding in different ambient densi- 
ties. Indeed, the data for SNRs in M 82 agree well with the ex- 
trapolation of the E-D relation derived for SNRs in other galax- 
ies to smaller sizes. Thus, the arguments we have exposed in this 
paper should also be applied to SNRs in M 82. 

Even though they have small sizes in the parsec range 
(Muxlow et al. 119941 ), they could be close to the end of their 
Sedov phase provided that the ambient densities are ~ 10 3 cm" 3 , 
with corresponding shock velocities of ~ 10 3 kms _1 com- 
pared to ~ 10 4 kms _1 as in the case of undecelerated expan- 
sion. These lower shock velocities infer characteristic times 
of ~ 10 3 yr, compatible with the average radio-emitting life- 
time found by Rronberg et al. (2000). The only exceptions 
are a few fast-evolving radio sources that are probably radio 
SNe, namely young objects still evolving in their circumstel- 
lar medium. Although the detailed physical conditions and pro- 
cesses at such high densities may differ from those at densities 
typical of other galaxies, we are confident that the simple esti- 
mate presented above is adequate to justify the large measured 
characteristic times within our framework. Finally, Fenech et al. 
(2008) claimed to have measured very high expansion veloci- 
ties, up to 10, 000 km s -1 : these measurements seem to conflict 
not only with our model but also with the above-mentioned pho- 



R. Bandiera & O. Petruk: A statistical approach to radio emission from shell-type SNRs 



9 



tometric measurements (Kronberg et al. |2000| l. Therefore, it is 
probably too early to discuss these latest measurements. 

6. Conclusions 

We have shown that studies of the statistical properties of SNR 
samples may provide insight into the physics of electron accel- 
eration and the time evolution of SNRs. We have proposed a 
new scenario, along the lines of previous work by Berkhuijsen 
( 1986 ), which interprets in a natural way the observed correla- 
tions between radio surface brightness S, size D, and ambient 
density n a in a sample of SNRs. The main parameter of SNR 
evolution that enters into these correlations is the time at which 
a SNR ceases to behave as a radio source, and we find that this 
endpoint is located close to the end of the Sedov phase. We oth- 
erwise find that the observed correlations mostly reflect that the 
sample consists of SNRs located in very different ambient condi- 
tions; while the evolution of individual SNRs plays a secondary 
role, and cannot be extracted by simply studying correlations be- 
tween pair of quantities. 

Within this framework, we present a new approach to ana- 
lyzing the statistical data, based on a 2-dimensional fit to 2 as a 
function of D and n . We show that the slope of 2(D) at constant 
n Q should represent more closely the true evolutionary track of 
an individual SNR than the well known "X-D relation", which 
is obtained without including information about n Q . 

In this work, we have used data published by Berkhuijsen 
( 1986 ). Although this data sample is rather limited, our method 
of analysis applied to these data is already capable of discrim- 
inating to some level between different theoretical models. For 
instance, models prescribing constant efficiencies for both mag- 
netic field (turbulent) amplification and electron acceleration 
(e.g., Berezhko & Volk 2004) are well outside the parameters 
region allowed by the data. On the other hand, models assuming 
a constant acceleration efficiency but a constant post-shock mag- 
netic field are marginally (about 2-<x) consistent with the data. 
For the sake of comparison, we have applied the same technique 
to a sample of SNRs in M 33. Although this sample could be af- 
fected by selection effects, it is completely independent from the 
sample of Berkhuijsen, and the parameters that we derive from 
the two samples are in close agreement, within the statistical er- 
rors. This may indicate again that our technique is robust andthe 
assumptions at the basis of our statistical analysis are correct. 
With the use of larger and more accurate data samples the con- 
fidence interval will become narrower, and the statistical results 
will be capable of constraining models more tightly. 

In nearby galaxies, deep and detailed surveys of SNRs have 
become available in both radio and X-rays. The number of SNRs 
in other galaxies has increased considerably, with reasonably 
good radio and X-ray flux measurements and in various cases 
also information about their angular size. Although SNRs in ex- 
ternal galaxies are more difficult objects to observe than Galactic 
SNRs, their distances are well known (since they correspond to 
the distance of the parent galaxy); in addition, selection effects 
can be modeled in a more accurate way for a sample of SNRs 
that are all at the same distance. Therefore, a natural extension 
of the present work is to elaborate our diagnostic tools along the 
guidelines shown in this paper, and to use them to investigate 
SNRs samples in various nearby galaxies. These results will be 
presented in a forthcoming paper. 
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